Feel++ : A Versatile High Performance Finite Element Embedded Library into C++

International audienceA Domain Specific Language for PDEs embedded in C++ providing a syntax very close to the mathematical language in order to describe the variational formulation: Supports generalized arbitrary order Galerkin methods (cG, dG) in 1D, 2D, 3D Supports simplex, hypercube and high order meshes Supports finite elements: Lagrange, Hermite, Nedelec, Raviart-Thomas Supports seamless parallel computing Supports seamless interpolation between grids/function spaces Supports symbolic computation thanks to GiNaC Supports large scale parallel linear and non-linear solvers (PETSc/SLEPc) Supports hybrid computing: MPI, multi-thread, GPU (HARTS) Supports in-situ visualization with ParaView</ul

Efficient solving strategies for incompressible Navier-Stokes equations for large scale simulations using the open source software Feel++

International audienceOver the past few decades, the computational fluid dynamics has evolved to become an important tool for the description of the complex multi-physics, multi-scale phenomena characterizing blood flow. Its reliability depends both on the verification of the numerical methods and on the validation of the mathematical models. The aim of the first part of the talk is to present a preconditioning framework for the linear system arising from the finite element discretizations and time advancing finite difference schemes of the 3D steady and unsteady Navier-Stokes equations. In particular, we are interested in preconditioners based on an algebraic factorization of the system's matrix which exploit its block structure, such as the Pressure Convection-Diffusion (PCD) preconditioner , the SIMPLE preconditioner or the LSC preconditioner, see [Elman et al. (2014)]. A comparison between the efficiency of these preconditioners is ascertain by testing them over the 3D backward facing step benchmark. The iteration counts using the PCD preconditioner are independent of mesh size and high order finite elements and mildly dependent on Reynolds numbers for steady flow problems which is not the case for the other preconditioners. In the second part of the talk we describe a framework for the solution of flow problems relevant to biomechanics strongly supported by the aforementioned solving strategies. We assess the efficiency of this framework through experimental data for fluid flow in a nozzle model with rigid boundaries, a device designed to reproduce acceleration, deceleration and recirculation, features commonly encountered in medical devices [Stewart et al. (2012)]. The flow rates were chosen to cover laminar (Re = 500), transient (Re = 2000) and turbulent (Re = 3500) regimes. The numerical results displayed during the presentation are obtained using the open-source library Feel++ (Finite Element method Embedded Language in C++, www.feelpp.org). Figure 1: Computational Fluid Dynamics FDA Benchmark at Re = 2000

Hemodynamic simulations in the cerebral venous network: A study on the influence of different modeling assumptions

International audienceBlood flow computations in complex geometries are of major interest in various cardio-vascular applications. However, deriving an appropriate computational model is still an open issue and a central question is how to incorporate and quantify uncertainties due to different modeling assumptions. The present work is intended as a first step in this direction, in the particular case of blood flow in the cerebral venous system. After a careful evaluation of the influence of the computational methodology, the study investigates the impact on the velocity field and the wall shear stress of three inflow boundary conditions, two strategies for treating the outflow boundary condition and two different viscosity models. The results demonstrate that the effect of setting the inflow boundary condition on the forces created by blood flow, is likely greater than for other modeling assumptions, the other important factor being the blood viscosity model, especially in wall shear stress computations. They suggest that improvements on the one hand on the mathematical and computational approach, and on the other hand on available data for their treatment are needed

Vers la simulation des écoulements sanguins

Contrairement aux liquides ordinaires, les fluides complexes comme le sang exhibent des comportements étranges qui dépendent essentiellement des structures sous-jacentes qui les composent. La simulation des écoulements sanguins continue de poser un formidable défi pour les modélisations théoriques et numériques dont l'intérêt est de développer des méthodes et des outils de simulation pour la communauté médicale. Nous proposons dans cette thèse une contribution à ce projet qui sera majoritairement centré sur les aspects numériques et informatiques. Nous nous sommes particulièrement intéressés à l'interaction entre le sang et la paroi vasculaire, qui joue un rôle important dans les grandes artères comme l'aorte. Nous nous sommes aussi investis dans la simulation du transport des cellules sanguines dans le sang. Pour la résolution des équations aux dérivées partielles décrivant nos modèles d'hémodynamique, nous avons choisi d'utiliser des méthodes numériques dont la précision pourra être accrue de manière arbitraire. Dans ce but, les principaux ingrédients qui ont été mis en oeuvre sont (i) la méthode des éléments finis basée sur des approximations de Galerkin d'ordre arbitraire en espace et géométrie, (ii) la méthode ALE pour la prise en compte de la mobilité des domaines pour des déplacements d'ordre arbitraire, (iii) les couplages implicites et semi-implicites pour l'interaction fluide-structure. Nous proposons également une nouvelle formulation de la méthode de la frontière élargie visant à modéliser le transport de particules déformables immergées dans un fluide. Nos simulations numériques se sont appuyées sur la librairie de calcul Feel++, spécialisée dans la résolution d'EDP. Outre l'implémentation des modèles physiques, nous y avons développé diverses fonctionnalités nécessaires à la mise en oeuvre de nos méthodes : interpolation, méthode de Galerkin non standard, méthode ALE, environnement pour l'interaction fluide-structure. De plus, de par la taille des géométries et la complexité des modèles mis en jeu, le passage au calcul parallèle a été indispensable pour pouvoir réaliser nos simulations. Ainsi, nous avons décrit le développement qui a été effectué dans cette librairie pour permettre le déploiement de nos programmes sur des architectures parallèles.Unlike ordinary liquids, complex fluids such as blood exhibit strange behavior mainly dependent underlying structures that compose them. Simulation of blood flow remains a formidable challenge for theoretical and numerical modeling whose interest is to develop methods and simulation tools for the medical community. We propose in this work a contribution to this project will be mainly focused on numerical and computational aspects. We are particularly interested in the interaction between the blood and the vascular wall, which plays an important role in the large arteries as the aorta. We also invested in transport simulation of blood cells in the blood. For solving partial differential equations describing our hemodynamic models we chose to use numerical methods whose accuracy can be increased arbitrarily. For this purpose, the main ingredients that have been used are textit {(i)} the finite element method based on Galerkin approximations of arbitrary order in space and geometry, (i) the ALE method for taking into account the mobility of areas for movements of arbitrary order, (ii) the implicit and semi-implicit coupling for fluid-structure interaction. We also propose a new formulation of the method of boundary extended to model the transport of deformable particles immersed in a fluid. Our numerical simulations were based on the library Feel++ specializing on the EDP resolution. In addition to implementing physical models, we have developed various features are necessary for the implementation our methods: interpolation, Galerkin non-standard method, ALE method, framework for fluid-structure interaction. In addition, the size and geometry of the complexity of the models involved, the transition to parallel computing was essential in order to achieve our simulations. Thus, we have described the development was carried out in the library to allow deployment of our programs on parallel architectures.SAVOIE-SCD - Bib.électronique (730659901) / SudocGRENOBLE1/INP-Bib.électronique (384210012) / SudocGRENOBLE2/3-Bib.électronique (384219901) / SudocSudocFranceF

From medical imaging to numerical simulations

International audienceIn the last 20 years there have been lots of progress in 3D medical imaging (such as Magnetic Resonance Imaging, MRI, and X-ray Computed Tomography, CT) and in particular in modalities to visualise vascular structures. The resulting images have been successfully used in various clinical applications, in particular for cerebrovascular pathologies (e.g., neurosurgery planning; stenoses, aneurysm or thrombosis quantification; arteriovenous malformation detection and follow-up, etc.). The complexity of the processing and analysis of these images (size, information vs noise, artifacts, etc) led to the development of imaging tools such as vessel filtering, segmentation and quantification. There is however, until now, no database of synthetic images and associated ground-truths (segmented data) available in cerebrovascular images contrary to morphological brain image analysis (e.g. brainweb).In the ANR Vivabrain project, we combine the skills of several communities: computer science, applied mathematics, biophysics, and medicine to remedy the aforementioned observation. In particular we focus on complex multi-disciplinary problems such as (i) the handling of inter-individual cerebrovascular variability, (ii) the generation of computational meshes, (iii) the simulation of blood flows in the complete cerebrovascular system 3D+time (3D+t) including calibration and validation and (iv) the accurate simulation of the physical processes involved in MRA acquisition sequences in order to finally obtain realistic virtual angiographic images

ANGIOTK : An Open Platform to reconstruct vessels from MRI images and simulate blood flows to ultimately provide Virtual Angiographies

National audienceFrom medical images to numerical simulations ANGIOTK The platform is still a prototype, but reached a level of maturity that requires very little human intervention. The overall objective is to make available, to the medical community, tools for computer aided modeling for conducting experiments in silico, otherwise difficult or impossible to carry on patients. A specific aim is to create an ecosystem around the platform composed of academic, research centers and enterprises and to provide services such as training, specific developments, consulting or deployment. The platform is operated by Cemosis and Kitware. An Open Platform to reconstruct vessels from MRI images and simulate blood flows to ultimately provide Virtual Angiographie

Towards Blood Flow Simulation

Contrairement aux liquides ordinaires, les fluides complexes comme le sang exhibent des comportements étranges qui dépendent essentiellement des structures sous-jacentes qui les composent. La simulation des écoulements sanguins continue de poser un formidable défi pour les modélisations théoriques et numériques dont l'intérêt est de développer des méthodes et des outils de simulation pour la communauté médicale. Nous proposons dans cette thèse une contribution à ce projet qui sera majoritairement centré sur les aspects numériques et informatiques. Nous nous sommes particulièrement intéressés à l'interaction entre le sang et la paroi vasculaire, qui joue un rôle important dans les grandes artères comme l'aorte. Nous nous sommes aussi investis dans la simulation du transport des cellules sanguines dans le sang. Pour la résolution des équations aux dérivées partielles décrivant nos modèles d'hémodynamique, nous avons choisi d'utiliser des méthodes numériques dont la précision pourra être accrue de manière arbitraire. Dans ce but, les principaux ingrédients qui ont été mis en oeuvre sont (i) la méthode des éléments finis basée sur des approximations de Galerkin d'ordre arbitraire en espace et géométrie, (ii) la méthode ALE pour la prise en compte de la mobilité des domaines pour des déplacements d'ordre arbitraire, (iii) les couplages implicites et semi-implicites pour l'interaction fluide-structure. Nous proposons également une nouvelle formulation de la méthode de la frontière élargie visant à modéliser le transport de particules déformables immergées dans un fluide. Nos simulations numériques se sont appuyées sur la librairie de calcul Feel++, spécialisée dans la résolution d'EDP. Outre l'implémentation des modèles physiques, nous y avons développé diverses fonctionnalités nécessaires à la mise en oeuvre de nos méthodes : interpolation, méthode de Galerkin non standard, méthode ALE, environnement pour l'interaction fluide-structure. De plus, de par la taille des géométries et la complexité des modèles mis en jeu, le passage au calcul parallèle a été indispensable pour pouvoir réaliser nos simulations. Ainsi, nous avons décrit le développement qui a été effectué dans cette librairie pour permettre le déploiement de nos programmes sur des architectures parallèles.Unlike ordinary liquids, complex fluids such as blood exhibit strange behavior mainly dependent underlying structures that compose them. Simulation of blood flow remains a formidable challenge for theoretical and numerical modeling whose interest is to develop methods and simulation tools for the medical community. We propose in this work a contribution to this project will be mainly focused on numerical and computational aspects. We are particularly interested in the interaction between the blood and the vascular wall, which plays an important role in the large arteries as the aorta. We also invested in transport simulation of blood cells in the blood. For solving partial differential equations describing our hemodynamic models we chose to use numerical methods whose accuracy can be increased arbitrarily. For this purpose, the main ingredients that have been used are textit {(i)} the finite element method based on Galerkin approximations of arbitrary order in space and geometry, (i) the ALE method for taking into account the mobility of areas for movements of arbitrary order, (ii) the implicit and semi-implicit coupling for fluid-structure interaction. We also propose a new formulation of the method of boundary extended to model the transport of deformable particles immersed in a fluid. Our numerical simulations were based on the library Feel++ specializing on the EDP resolution. In addition to implementing physical models, we have developed various features are necessary for the implementation our methods: interpolation, Galerkin non-standard method, ALE method, framework for fluid-structure interaction. In addition, the size and geometry of the complexity of the models involved, the transition to parallel computing was essential in order to achieve our simulations. Thus, we have described the development was carried out in the library to allow deployment of our programs on parallel architectures

Vers la simulation des écoulements sanguins

Unlike ordinary liquids, complex fluids such as blood exhibit strange behavior mainly dependent underlying structures that compose them. Simulation of blood flow remains a formidable challenge for theoretical and numerical modeling whose interest is to develop methods and simulation tools for the medical community. We propose in this work a contribution to this project will be mainly focused on numerical and computational aspects. We are particularly interested in the interaction between the blood and the vascular wall, which plays an important role in the large arteries as the aorta. We also invested in transport simulation of blood cells in the blood. For solving partial differential equations describing our hemodynamic models we chose to use numerical methods whose accuracy can be increased arbitrarily. For this purpose, the main ingredients that have been used are textit {(i)} the finite element method based on Galerkin approximations of arbitrary order in space and geometry, (i) the ALE method for taking into account the mobility of areas for movements of arbitrary order, (ii) the implicit and semi-implicit coupling for fluid-structure interaction. We also propose a new formulation of the method of boundary extended to model the transport of deformable particles immersed in a fluid. Our numerical simulations were based on the library Feel++ specializing on the EDP resolution. In addition to implementing physical models, we have developed various features are necessary for the implementation our methods: interpolation, Galerkin non-standard method, ALE method, framework for fluid-structure interaction. In addition, the size and geometry of the complexity of the models involved, the transition to parallel computing was essential in order to achieve our simulations. Thus, we have described the development was carried out in the library to allow deployment of our programs on parallel architectures.Contrairement aux liquides ordinaires, les fluides complexes comme le sang exhibent des comportements étranges qui dépendent essentiellement des structures sous-jacentes qui les composent. La simulation des écoulements sanguins continue de poser un formidable défi pour les modélisations théoriques et numériques dont l'intérêt est de développer des méthodes et des outils de simulation pour la communauté médicale. Nous proposons dans cette thèse une contribution à ce projet qui sera majoritairement centré sur les aspects numériques et informatiques. Nous nous sommes particulièrement intéressés à l'interaction entre le sang et la paroi vasculaire, qui joue un rôle important dans les grandes artères comme l'aorte. Nous nous sommes aussi investis dans la simulation du transport des cellules sanguines dans le sang. Pour la résolution des équations aux dérivées partielles décrivant nos modèles d'hémodynamique, nous avons choisi d'utiliser des méthodes numériques dont la précision pourra être accrue de manière arbitraire. Dans ce but, les principaux ingrédients qui ont été mis en oeuvre sont (i) la méthode des éléments finis basée sur des approximations de Galerkin d'ordre arbitraire en espace et géométrie, (ii) la méthode ALE pour la prise en compte de la mobilité des domaines pour des déplacements d'ordre arbitraire, (iii) les couplages implicites et semi-implicites pour l'interaction fluide-structure. Nous proposons également une nouvelle formulation de la méthode de la frontière élargie visant à modéliser le transport de particules déformables immergées dans un fluide. Nos simulations numériques se sont appuyées sur la librairie de calcul Feel++, spécialisée dans la résolution d'EDP. Outre l'implémentation des modèles physiques, nous y avons développé diverses fonctionnalités nécessaires à la mise en oeuvre de nos méthodes : interpolation, méthode de Galerkin non standard, méthode ALE, environnement pour l'interaction fluide-structure. De plus, de par la taille des géométries et la complexité des modèles mis en jeu, le passage au calcul parallèle a été indispensable pour pouvoir réaliser nos simulations. Ainsi, nous avons décrit le développement qui a été effectué dans cette librairie pour permettre le déploiement de nos programmes sur des architectures parallèles